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Multi-nucleon correlations in Deep Inelastic Scattering at large Bjorken xb 
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Using realistic nuclear many-body wave functions in the collinear factorization framework, we 
evaluate nuclear Deep Inelastic Scattering at large Bjorken xb, where correlations are needed to 
exceed the kinematic limits of a single free nucleon. We find that recent experimental data on cross- 
section ratios from the CLAS collaboration exhibit features that are not reproduced by modern 2- 
and 3-nucleon correlations in the nuclear wave function. New physical mechanisms are needed to 
explain data in the 2 < xb < 3 range. 
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Introduction. The CLAS collaboration has recently 
published data on ratios of cross-sections for scattering 
a high-energy electron on a nucleus A to 3 He at large 
xb < 3 [H, Q- These ratios show an interesting structure 
with 2 plateaus at xb > 1.5 and xb > 2.25, of height 
increasing with A. Earlier SLAC data on A/ 2 H ratios 
3) and A/ 4 He ratios [3| exhibit features similar to those 
seen by CLAS. The emergence of plateaus is generally 
interpreted as due to 2- and 3-nucleon short-range corre- 
lations in nuclei. However, the physics underlying these 
correlations has remained elusive. Short-distance multi- 
quark [f| or multi-nucleon clustering @, H, 0] has been 
favored in several models. See 7, 8] for a review. 

In this paper we demonstrate the limited extent to 
which the observed plateaus can be explained in terms 
of realistic inter-nucleon potentials including local short- 
distance repulsion 0, [lfj |. which can push a nucleon 
momentum to large values. The resulting 2-nucleon 
(NN) correlations produce a nearly nucleus-independent 
shape of single nucleon momentum I distributions in the 
2 I ^$ 4 fm _1 range 0, EH, which may explain the 
first plateau. At higher momenta, 3-nucleon (NNN) cor- 
relations become dominant Q. We investigate whether 
the NNN correlations may explain the second plateau. 

Our main theoretical inputs are realistic nuclear wave 
functions [l(J, and the collinear factorization frame- 
work to describe the scattering of the lepton on a bound 
nucleon in terms of electron-parton scatterings. The ob- 
tained formulas are rather general and the few approxi- 
mations are made explicit. The formalism is easily gen- 
eralizable to include other physical mechanisms. 

We provide a minimalist approach that includes only 
the effect of "Fermi motion" , or single-nucleon momen- 
tum distributions, with exact treatment of kinematics 
at the parton, nucleon and nuclear level. It provides a 
baseline, including only nucleon degrees of freedom and 
the best available nuclear wave- functions. We find that, 
while the first plateau may be explained in our model, 
the second plateau in the 2 < xb < 3 region is not de- 
scribed by inclusion of NNN correlations. Thus, single- 



nucleon degrees of freedom are inadequate, and more ex- 
otic physical mechanisms such as those mentioned above 
are needed to explain experimental data. 

Collinear factorization. Under the one-photon ap- 
proximation, the cross section of nuclear deep inelastic 
scattering is determined by the nuclear hadronic tensor, 

, A ,q) = ±- [ (Pze-^frAlj^WmpA) , (1) 



47T 



with nuclear momentum Pa = ApA, virtual photon 
momentum q, and corresponding current j. At large 
Bjorken xb — Q 2 /2pa ■ q with Q 2 = —q 2 , the virtual 
photon is very localized in space-time. In the impulse 
approximation, the virtual photon interacts with a par- 
ton of momentum k belonging to a nucleon of momen- 
tum p, bound in a nucleus of momentum Pa and mass 
Ma = Am, as sketched in Fig. [1] We picture the nu- 
cleus as a system of A bound nucleons whose momen- 
tum distributions ( "Fermi motion" ) are given by nuclear 
many-body wave functions. If we assume the Hilbert 
space of the nucleus to be the direct product of A single- 
nucleon Hilbert spaces, and the current operator to act 
on a single-nucleon subspace, 

W£ v (p A ,q)* [ d^zl- [ dn A e^ z (p\3^(z)f(0)\p) 
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= / dfi A W^ u (p,q) 



(2) 



with the hadronic tensor of a bound nucleon, and 
relativistic Fermi smearing measure, with m 2 = p 2 , 

d ^ A = ^W^ PAip \o=V^ ■ (3) 

The single nucleon 4-momcntum distribution (ND) is 



Pa(p) 



n 

i=2,A 



l(P,P2:--->Pa)| 5 



(4) 
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FIG. 1: Forward Compton amplitude for nuclear DIS. 



with 4>a the nuclear wave-function. Gauge invariance and 
Eq. ([2|) imply q^W^ — q^W^ — 0, so we can define the 
nuclear and nucleon structure functions as follows: 



W^(x B ,Q 2 ) = -r v F 1A {x B ,Q A ) 



_ PaPa 

PA ■ q 

pup" 



F 2A (x B ,Q 2 ) 



q^q v /q 2 , and 



WT(^.0 ) = -aTFiixN, Q z ) + ^F 2 (x N , Q z ) (5) 

where xn — Q 2 /2p ■ q and g^ — - 
a M = a M — (a • q/q 2 )q^ for any 4-vector a. 

To evaluate the hadronic tensor or the structure func- 
tions in terms of the QCD collinear factorization ap- 
proach [13], we need to define suitable light-cone "+" 
and "— " directions. We have 2 choices: using the (q,p A ) 
plane (A-frame) or using the (q,p) plane (N- frame). Our 
goal is to express the nuclear structure functions in terms 
of nucleon structure functions. Hence, as we will discuss, 
the optimal choice is the N-frame, 



p n H —n' 

2p+ 

-Uup+n^ + 



(6) 



p^ A = uip + n fl 



2uip 



where n = (1/V2, 0±, l/Vty, n = [1/V2, 0±, -1/V2), 
and = (do ± o,s) In Eq. 



^ = Pa/p + is 



the inverse nucleon fractional plus-momentum with re- 
spect to the nucleus, = m 2 + p\ ± , and £a — 

2xb/(1 + \ 1 + 4:X 2 B m 2 j_/ Q 2 ) is the nuclear Nachtmann 
variable [13[ in the N-frame. 

Following QCD collinear factorization, we expand the 
parton momentum k, entering the top blob in Fig. [TJ 
around its positive light-cone component, fc M = xp + n tl + 
0(k — xp + n) with parton momentum fraction, x = 
k + /p + , and factorize the nucleon hadronic tensor as 

W£»(x N ,Q*) = J2 J ^K? (x,Q 2 ) Vf/N (x,Q 2 ) 



0(1/Q 2 ) 



(7) 



where x = ^a^/(1— (Caw) 2 m 2 /Q 2 ) and x — Q 2 /2k-q = 
{£,au)/x in the N-frame. In Eq. ||7J), <fif/N is the leading 
twist PDF of a bound off-shell nucleon for a parton of 
flavor / = g, q, q, and Ti^" is the partonic tensor for 
an on-shell parton of momentum fc M = xp + n^ with all 
perturbative collinear divergences along the parton mo- 
mentum k absorbed into the PDFs. With no surprise, 
Eq. (f7|) has the same factorized form as that of a free 



nucleon [lj] because the factorization of short-distance 



partonic dynamics is insensitive to the details of long- 
distance hadron physics. From the parton level gauge 
invariance, q^Ti.'j 1 ' = 0, we can decompose the partonic 
tensor as 

Hf(x, Q 2 ) = -gT h){x, Q 2 ) + ^ h%x, Q 2 ) , (8) 

where the scalar functions can be computed order-by- 
order in powers of a s from the factorized expression (J7|) 
with the nucleon state N replaced by a parton state of 
flavor /, regardless of the nucleon state. 

Combining Eqs. ©, ©, (0 and §8§ we derive nuclear 
structure functions in terms of leading twist massless 
(to 2 = 0) nucleon structure functions, with i = 1, 2 



F 2A {x B ,Q 2 ) 



XB 



l + S A 



3(1 + 5 U ) 2 



(1 + <Sa)(1+*») 



F 1A (x Bl Q 2 ) = j d^F[ 0) {Uu,Q 2 ) 



(1 



f 2 (0) (£aw,Q 2 ), 



2U" 



(1 + 8 a )(1 + S„) 



f!?\uu,Q 2 ) 

4£a^ 



6(1 - U") 



(9) 
(10) 



where 5a — 4,x 2 B m 2 /Q 2 , 5 n — Ax 2 N m 2 /Q 2 , 5 U = 
4x B x N M 2 /Q 2 , and M 2 = f [m 2 + m\/uj 2 } . 

We would like to emphasize that only in the N-frame 
we can easily generalize the free-nucleon PDF to the PDF 
of a bound off-shell nucleon because they both have the 
same operator definition. For example, other than the 



state \p), the definition of quark PDF at leading order, 

is the same for both bound and free nucleon This 
would not be true in the A-frame. 
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Large- xb correlations. To apply the general result, 
Eqs. (|9j)- (fT0l) . to describe CLAS data we need to choose 
of bound nucleons and the measure dp a- First, 
we assume the nucleons to be on their mass-shell with 



??? 



2 _ ^2 



m 2 . Thus, we may identify with the massless 
free nucleon structure functions. They can be computed 
using PDFs from QCD global fits which do not already 
correct for the target's mass. Corrections for off-shell nu- 
cleons were discussed in Ref. [l6| in a related, but some- 
what different, formalism. Next, we use non-relativistic 
computations of nucleon distribution functions 

[lEol (we 

are not aware of any relativistic model which includes NN 
and NNN correlations). In summary we approximate 



p A (p) * (2n) 4 2p 5(p 2 -m*) pj(p) 



(12) 



Defining p A = p + I, using the translation invariance of 
the nucleon distribution, and choosing the struck nucleon 
rest frame, p = (m, 0), we obtain dpA = d 3 lp* A (l) , where 

p*a(i) = JOIUa d%)\Mh h,... Ja)\ 2 s^ (i + Ef =2 U) 

can be identified with the nucleon distribution computed 
in the nucleus rest frame. We consider state-of-the- 
art non-relativistic nucleon distribution functions p* A ob- 
tained in a Variational Monte Carlo (VMC) computation 
which include NN and NNN correlations [9( , compared to 
the parametrization of p* A discussed by Ciofi degli Atti 
and Simula (CS) [1(1, which only considers NN correla- 
tions, see Fig. [2] The relative momentum I = (l±,ls) 
further enters Eq. © through to = (I3 + y/f-i- rri\_)/m 
and rn\ = rn 2 + 1\ . The per- nucleon nuclear deep inelas- 
tic (DIS) cross section reads 



da A 



4na 2 \ 1 



dQ 2 dxi 



-y 2 F 1A (x B ) 



Q 4 U 

9 

m 2 , 
1-V-qS*bV 



F 2 a(xb) 



XB 



(13) 



where y — Pa - q/PA • Pi = Q 2 /(2Ei a bmxB), Pi is the 
lepton initial 4-momentum and Ei a b its energy in the 
target rest frame. The nuclear structure functions FiA 
are given in Eqs. (j9|)- ([T0|l . For the following plots, we 



compute i^ (0) at leading order using CTEQ5L PDF [17 1. 



In Fig. [3] left, we compare the A/ 3 He ratios from 
Eq.(H2]) with CLAS data taken at Q 2 ~ 2 GeV 2 [J. After 
an initial rise, the computed ratio becomes softer start- 
ing at xb > 1-5. This is due to the onset of NN correla- 
tions which develop a tail in the nucleon distribution at 
|^| > 2 fm -1 , see Fig. [5J Our curves using VMC distri- 
butions (solid lines) can describe reasonably well CLAS 
data for light nuclei. However, they do not describe the 
onset of the second plateau, otherwise visible in the ND 
ratios of Fig. [2] and due to the onset of NNN correla- 
tions. Indeed, in the Fermi smearing of the structure 
functions ([9j)- (fT0|) . the nucleon momentum I is sampled 
with a rather large standard deviation of about 1 fm -1 , 
which smears out such structure. This second plateau 




FIG. 2: Top: VMC nucleon momentum distributions [jj. 
Bottom: A/ 3 He ratios (thick: VMC; thin: CS [lC 



is then suggestive of new mechanisms such as 6- and 
9-quark clusters [§:.], multi- nucleon clusters [3, 0, 0], or 
dynamical short-range nucleon correlations [l8| , effective 
when nucleons approach each other at a distance smaller 
than a nucleon radius. In light nuclei such as 3 He or 4 He, 
all nucleons exist on the nuclear surface, which reduces 
the probability of nucleon overlap. In heavier nuclei, an 
increasing fraction of nucleons exist in the nuclear inte- 
rior where the overlap probability is maximized, so that 
we expect the difference between our baseline and exper- 
imental data to increase with A at x > 2.25. The com- 
putation with the CS distributions (dashed lines) fails at 
xb > 1-5 for nuclei heavier than 4 He, perhaps due to the 
lack of NNN correlations. 

In Fig. [3] right, we show the effect of Fermi motion on 
A/ 3 He ratios at Q 2 = 5 GeV 2 . The onset of correlation 
effects has moved down to xb ~ 1 and the plateau's 
slope has increased. At still higher values of Q 2 we found 
that the onset of correlation effects does not move much 
further, but the slope of the plateau continues to increase. 

We estimated the sensitivity of the A/B ratios to dif- 
ferent choices of PDF, which is at most 5%. Parton 
recombination effects, which are important at small-xs 
0, are negligible in our case. 

The simple model of the nucleus we have considered 
can be improved in many ways. First, we used a non- 
relativistic computation of pa to describe the large xb re- 
gion, to which nucleons of relative momentum I > 2 fm -1 
contribute: one may expect relativistic corrections to be- 
come important and change the hard tails of the ND. Sec- 
ond, the nucleon hadronic tensor can have a more general 
Lorentz decomposition than in Eq. ([5]), as discussed in 
[l6| where our assumption leading to Eq. @ was relaxed. 
Finally, we have made no attempt to describe the EMC 
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FIG. 3: Left: Comparison of cross-section ratios computed at Q 2 = 2 GeV 2 and Ei a b = 4.5 GeV to CLAS data taken at 1.6 
GeV 2 < Q 2 < 2.4 GeV 2 [l[. Solid line computed with VMC nucleon distributions dashed line with the CS parametrization 
P1J. Right: cross-section ratios at Q 2 = 5 GeV 2 and -E; a f> = 9 GeV, accessible at the 12 GeV upgrade at JLAB. 



effect in the 0.2 < xb < 1 region, see 2(J 21 1 for a review. 
Proposed underlying mechanisms include off-shell correc- 
tions to nucleon structure functions, nucleon binding and 
removal energy corrections, and Q 2 -rescaling. Their im- 
pact at large- a; b needs to be estimated, but we do not ex- 
pect major changes to our results since the cross-section 
ratios at large xb are mainly dominated by the relative 
strength of the ND tails. Quasi-elastic lepton- nucleon 
scattering, which causes the dip in experimental data at 
xb ~ 1 can be straightforwardly implemented, as well 
as contributions from the nuclear pionic cloud. However, 
such improvements are not expected to make significant 
contributions at large xb- 

Finally, at the lower Q 2 we considered, a few theoreti- 
cal issues arise. First, especially at Q 2 as 2 GeV 2 , one can 
question the partonic picture of DIS for computing F^a- 
When computing the ratio, however, this should not be 
an issue thanks to parton hadron duality. Second, one 
would expect higher-twist corrections to become sizable. 
They would not completely cancel in the ratio, because 
they are expected to be proportional to A 1 / 3 . We will 
consider these issues in a future work, but we expect 
these improvements to provide cross section ratios with 
smooth corrections that will not likely induce the second 
plateau structure seen in the CLAS data. 
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